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Electron scattering in the monolayer graphene with short-range impurities modelled by the annu- 
lar well with a band-asymmetric potential has been considered. Band-asymmetry of the potential 
resulted in the mass (gap) perturbation in the Dirac equation. Exact explicit formulae for the scat- 
tering matrix have been derived. The results are presented in terms of the scattering phases and in 
the geometrical form of a relation between some 2-vectors. The characteristic equation is obtained. 
It has a form of the orthogonality condition. An approximate calculation of observables in terms 
the scattering theory results is outlined. 



I. INTRODUCTION 



We consider in this work the electrons scattering in the 2-1-1 Dirac equation model of the monolayer graphene due 
to the short-range perturbations. Short-range potential impurities in graphene were considered in works [1], 
In our work [3| , a new model of the short-range impurities in graphene was considered taking into account the obvious 
fact that the Kohn-Luttinger matrix elements of the short-range perturbation calculated on the upper and lower band 
wave functions are not equal in a general case. This means that the perturbation must be generically described by 
a Hermitian matrix. The diagonal matrix case corresponding to a presence of the potential and mass perturbation 
was considered. Annular geometry of the perturbation allows us to exclude a perturbation at small distances, which 
are unphysical in a crystal. Delta function model of the short-range diagonal matrix perturbation was considered 
in our works [1], We study here the potential and mass well model. We consider scattering of electrons by this 
potential and obtain exact explicit formulae for the scattering data. We also discuss application of these formulae for 
calculation of observables. 



II. BASIC EQUATIONS 

The Dirac equation for electronic states in graphene for our model of matrix annular well described above reads 



\^-ihvF a^d^ -P{m + 5m [r)) j ij = {E-V (r)) V, (1) 

where vp is the Fermi velocity of the band electrons, /3 are the Dirac matrices /3 = (T3, ai = cri, q;2 = i<^2i are 
the Pauli matrices, 2mfF^ = Eg is the electronic bandgap, il) (r) is the two-component spinor. The electronic gap can 
appear in the graphene monatomic film lying on the substrate because of the sublattices mutual shift fTj. The spinor 
structure takes into account the two-sublattice configuration of graphene. 5m (r) and V{r) are the local perturbations 
of the mass (gap) and the chemical potential: 

V (r) = -aA (r) , 5m (r) = -feA (r) , (2) 

where A (r) is determined as follows 

0, < r < ri (region I), 
A (r) = 1, Ti <r <r2 (region II), . (3) 
0, r >r2 (region III). 
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Let us present the two-component spinor in the form 

exp(-zM) /7,Mexp[z(j-l/2)^] 



9j W exp[i(j + 1/2) (p] 



(4) 



where j is the pseudospin quantmxi number; j ~ ±1/2, ±3/2, . . .. In opposite to the relativistic theory, this quantum 
number has nothing to do with the real spin and indicates a degeneracy in the biconic Dirac point. The upper fj (r) 
and lower gj (r) components of the spinor satisfy the equations set 



(5) 



^ + -lf,-{E + m)g, = ia-b)A{r)g,. 



These equations have a symmetry: 



-J, a 



(6) 



(7) 



III. SCATTERING MATRIX AND CHARACTERISTIC EQUATION 

When r'e [ri, (in regions I and III), we obtain from eq. (O and eq. 



p2 2 i (i - 1) 

E — m 



= 0. 



This equation is related to the Bessel one. Its solution in the region I reads: 



where k = ^/E"^ — m^, J„ (x) is the Bessel function. 
Let us introduce the function ip^ (r) : 

We obtain from ([9]) and eq. ([T0|) for the region < ?' < ri: 



E + m Jj-i/2 {^.r) 
E - m J,+i/2 [kt) 



We have from ([3]), ([S]), and © within the region II: 

'^^^ +lg,-{E-m)f, = {a + b)f,, 



dr 



(8) 

(9) 
(10) 

(11) 
(12) 

(13) 



^ + lf^-(E + m)f,^{a~b)g,. 



These equations can be re-written in the form 



f + '-f,-{E + fh]g,^0^ 



(14) 

(15) 
(16) 
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where E = E + a, m = m — b. Thus we have the following second-order equation in the region ri < r < r2 '■ 



,A = 0. 



Then the function ipj (r) (fTTj) can be written in the region II as follows: 

//, . E + m (fir) + (kt) 

^''^'^ E-m Jj+i/2 (kt) + (Kr) 

where — E"^ — rh. Similarly we obtain for the 



region III: 



(17) 



(18) 



(19) 



Continuity of the spinor components leads to matching conditions for the functions Lpj (r) at the boundaries between 
regions I, II and III. We obtain in result the following expressions for the coefficients Cj and Sj : 



'~"Nj_i/2 (Kri) Jj+i/2 {uri) - y/ 15^^^+1/2 (firi) Jj_i/2 (^ri) 



E+m 
E-fh 



F, 



(2) 



(20) 



(21) 



where 



\f^^H'f\/^{Kr2) [Jj+i/2 {nr2) + CjiVj+1/2 {nr2)\ 



IE 


■f m 




— m 


Ie- 


f m 



E — m 



(22) 



The constant Cj is determined by the formula (PU)) . Here a takes values 0, 1. The constant Sj is a phase factor of 
the outgoing wave, i. e. the S-matrix element in the angular momentum representation. Since Hn"^ (z) = Hn^* (z) 
for real z, the scattering matrix is unitary everywhere on the continuous spectrum. Equations (pill) . (pij) and ([^^ 
solve the electron scattering problem for the given potential. The denominator of Sj{E) is just the left-hand side of 
the characteristic equation for the bound and resonance electronic states: 



f(i) = 



(23) 



or 



^±^H^'] {nr2) [jj+1/2 {nr2) + C,iV,+i/2 {Kr2)] - \I^^Hf\^^{Kr2) [Jj_i/2 («r2) + C,N^_^,2 («:r2)] = 0. 

(24) 

Excluding the constant Cj from eq. (I^TI) and eq. (IMl) , we obtain the explicit formula for Sj and the characteristic 
equation: 

^/|^J,-+i/2 (^ri) J,-„i/2 (^ri) - v/|^J,--i/2 (^ri) J,+i/2 (^^ri) 
^i±SiV,+i/2 (^ri) J,_i/2 (^ri) - yji|iV,._i/, (^5.2) J,+i/2 (^n) 



1/2 (kJ-i) Jj-1/2 (^ri) - y|±^iVj_i/2 (Kr2) Jj+1/2 (Kri) 

_ v/|^'^.-+l/2 (^^2) gji\/2(^^2) - TM-^^- (^"^) ^j+\/2(^^2) 
y|iiiV, + i/2 (^r2) i?ji\/2(«r2) - y|±i^^-V2 («'^2) i^]+\/2(-'^2) ^ 



(25) 



When ri = 0, r2 = rp, we have a case of the simple round well; the characteristic equation takes the form: 



' E + m 

E ~ m 



Jj-i/2 (kto) ,2(Kro) 



E + m 
E — m 



Jj+i/2 {htq) H^^\,^iKro) 



(26) 



A geometric interpretation of these expressions for the S-matrix and characteristic equation can be done here. Let 
us introduce the following column vectors: 




(27) 



\ 



\ 



E+m j^(a 



(28) 



J-l/2 



{nr2) j 



and the matrix 



(r) 



Then ea.(P(7| can be written as follows 



Nj-i/2(iir) Nj+ii2{nr) 



(29) 



(30) 



where the vector Jj 




is determined as follows 



J,{r,)^D,{r,)3,{r^). 
r(") 



(31) 



Similarly, introducing the transformed vector H^"^ (r2) — D (r2) h^") (r2) , we can wtite the S-matrix ([2T|) . (|22|) in a 
simpler form: 



r(l) 



This formula can be written in another form introducing the components transposition operator P — ai : 

(5,(ri)J,,i(ri),ai5, (r2)hf (r2)) 



(32) 



Sj = - 



(5, (ri)J,,i(ri),ai5,(r2)hW (r2)) ' 
where (a, b) is a scalar product of the vectors a and b. The formula (I33p can be re-written in the form: 



(33) 



i^,J,(ri),hf) (r2) 



(if,J, (ri),hW (r2))' 



where the matrix Ki is determined as follows 



The characteristic equation now reads 



(i?,J, (ri),hW (r2) 



= 



(34) 

(35) 
(36) 
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IV. S- AND T-MATRIX PROPERTIES. POSSIBLE APPLICATIONS 

Using the relations Hn^ (z) = J„ + iNn, Hn^ ~ J„ — iNn, we can write S-matrix in the form: 

^ A,iE)+zB,iE) ^ B,{E)+iA,{E) 
^ A,{E)-zB,{E) B,{E)-iA,{Ey ^ ' 

and, therefore, it can be presented in the standard form [loj 

(E) = exp [i26j {E)] , (38) 

where the scattering phase is given by the expression 

<5,(i?) = arctan|l|||. (39) 

Formulae (j37l) . psp show once more that the scattering matrix Sj (E) is unitary on the continuum spectrum. The 
functions Aj (E) and Bj (E) are determined as follows: 



j E ~\- m ^ ^ E A- in ^ ^ 

^ = \ {'^^2) [Jj+1/2 (^ra) + (Krs)] - \ ^ -Jj+1/2 ('"^2) [Jj-1/2 (^ra) + (/tra)] 

V£/ — m \ E — m 



E — m 

(40) 



B = y ~'^T ^j+i/2 {^^^2) [Jj-1/2 (^ra) + CjiVj_i/2 (^ra)] - J ^^^ Nj_i/2 (wra) [Jj+i/a (^ra) + QiVj-+i/a (wra)] 
\ E — m \ — m 

(41) 

where the constant Cj is given by eq. ^0^. It is seen from ((201), (gOl), and (gl]) that all Sj (E) (j==±l/2, ±3/2 . . .) 
vanish , when a and b tend to zero, i. e. in the absence of a perturbation. 

It is easy to show that the phase is proportional to nro in the long-wave limit as it is necessary 

0, i (here 

''0 ~ fij ^2) (see also Q). The scattering amplitude / {9) and transport cross-section Yitr can be expressed in terms 
oiSj [E) as follow 0: 

f{e)^-^ [S,{E)-l]eMi{j- 1/2)9], (42) 

*^"'^'*i=±l/2,±3/2,... 



T.tr = 2/K ^ sin^ ((5j+i - (5j) (43) 

J = ±l/2,±3/2,.. 

In the vicinity of the resonance state energy, the Breit-Wigner form of the phase is valid [l3| : 

5, « 5f + arctan ^ , (44) 

2 [Ef e) 

where e'^^ and Tj are respectively the position and width of the resonance level, S'^p is the slowly-varying potential 
scattering phase. The presented above formulae can be used in order to calculate the Boltzmann conductivity JJj]: 

-^'^tr, (45) 

where the transport relaxation time equals 



l/rt,. = NiVF^tr- (46) 

Here iV^ is the impurities areal density, Ep = vpKp. The above equations transform the scattering data into the 
correspondent dependence of the Boltzmann conductivity. Thus characteristic features of the scattering data determine 
a behaviour of the electric conductivity. Proper numeric calculations will be presented elsewhere. 



6 



On the other hand, the derived here formulae for scattering data can be used to obtain an approximate density of 
states and other observables in fower order to the impurities density. The exact expHcit formula for one-impurity S- 
matrix obtained here, which allows us to calculate the scattering amplitudes, bound and resonance states, determines 
also the on-shell one- impurity T- matrix [l3j : 

T;- {E) = {l/i) [S, [E) - 1] (47) 

Corresponding off-shell T-matrix T°^^ (k, k',i?) (here k generically does not lie on the mass shell, i. e. \s? + m? can 
be not equal to E'^) can be written in the form [T2| : 

T°fi (k, k', E) - (vI/o(k)f {E) vPo (k')) . (48) 

where the transition operator T°f^ (E) is determined by the operator equation p^ : 

fo// ^ [/_ [/Gof°//. (49) 

Here Go = Hq^ is the operator inverse to the free Hamiltonian. All components of the product in eq. (|48l) are 
determined at different energies: E, E (k) and ECk'). Let us write the Lippmann-Schwinger equation (|49|) for off- 
shell one-impurity T-matrix in the matrix form: 

t" (k, k',E)^U{k-k')- f d\V (k - q) Y'X'V.'^-j " (q- k', i?) , (50) 

J q -f — E'' + lO 

where the Fourier transform of the potential is determined as follows 

U (k-k') = / (fr exp [-i (k - k') r] U{r) = 27r V e,„ cos m0 [ drrJ,n {kr) U (r) J„ (fc'r) , (51) 

1 if m = 0, , . 

2 ifmT^O. ' ' 

7p — fia^, Jo — /3, 9 is the angle between the vectors k and k'. 

Now we consider a question: how the obtained above formalae for scattering data can be used for calculation of 
observables in the monolayer graphene. It is seen from ([5T|) that in the case of the short-range perturbation, U (q) is 
a slowly- varying function. Then equation (jSOp shows that the off-shell T-matrix is a slowly- varying function of k and 
k' as well. At the same time, it can be a sharp function of energy E in the vicinity of a resonance. Therefore, in the 
case of a narrow resonance, the most important information on the off-shell transition matrix is given by the on-shell 
matrix, which is determined by the S-matrix (see (j47|) ). 

t"' (k, k', E) ~ T " (k, k', E)^2iY, ej-i/2 cos [(j - 1/2) 9] {S, {E) - 1} . (53) 

i 

Using the derived above explicit formula for S-matrix ([M)) (or (HI])), an approximate formula for the off-shell one- 
impurity T-matrix in the momentum representation can be obtained. Substituting the partial-wave projection of the 
S-matrix psp into equation (j53p . we can obtain an expression for the T-matrix, which can be used for approximate 
calculation of such observables as the denstity of states, for instance. The scattering phases 5j (E) are completely 
determined by the formulae (|40| , (|4T|) and (p9| obtained in this paper . The electronic density of states and other 
observables for systems with non-overlapping impurity potentials can be calculated on the base of Lloyd's formula 
(see, for instance, fl^). The simplest approximation can be obtained in the lower-order approximation to the mass 
operator M (k, E) in the impurities density n using the mean T-matrix method for a random distribution of impurities 

m 

M (k,!!,) ^n{k\T°^^ {E)\k) (54) 

Thus, having the above formulae for the scattering data, we can calculate the averaged over impurities distribution 
Green functions, density of states and the optical absorption coefficient. Results of such calculations will be published 
elsewhere. 
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V. CONCLUSION 



We considered the electron scattering problem in the monolayer graphene with short-range impurities. Charac- 
teristic for the electronic two-band theory band asymmetry of the potential is equivalently described by the scalar 
potential and the mass (gap) local perturbation. The crystal perturbation by a single impurity is modelled by the 
annular well. Exact explicit formulae for a single-impurity S-matrix and other scattering data have been obtained 
and analyzed for the short-range perturbation modelled with a use of the annular well with the band-asymmetric 
potential. The characteistic equation for bound and resonance states is derived for the case. This equation was shown 
to have a form of 2-vectors orthogonality condition. Possible application of these results to description of kinetical 
and optical properties of graphene is discussed. A procedure of approximate calculation of observables based on a 
substiution of the off-shell T-matrix by the one-impurity on-shell T-matrix, which can be expressed in terms of an 
exact S-matrix in the vicinity of a sharp resonance is suggested. The obtained in this paper explicit formulae for 
S-matrix can be used in this procedure. As an elementary example of this approach, we considered an approximate 
calculation of the mass operator in graphene with non-overlapping defects. 
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